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DATA.  REDUCTION  FOR  STABLE  AURORAL  RED  ARCS  OBSERVED  AT 
RAPID  CITY,  SOUTH  DAKOTA 


-x- 
J.  E.  Cruz,  Robert  Davies  ,  L.  K.  Droppleman,  E.  Marovich, 

+ 

L.  R.  Me gill,  M.  H.  Rees  ,  Linda  Reisbeck,  and  F.  E.  Roach 


All  sky  photometric  surveys  in  the  6300  A  [Ol] 
radiation  obtained  during  the  International  Geo- 
physical Year  at  Rapid  City,  South  Dakota  have  been 
reduced.   The  analysis  is  described  and  the  results 
for  the  night  of  September  16/17,  1958  are  printed  as 
isophote  maps.   These  maps  clearly  show  the  pres- 
ence of  an  arc  and  can  be,  therefore,  a  useful  tool 
for  further  investigations  of  the  stable  auroral  red 
arc. 

INTRODUCTION 
During  the  International  Geophysical  Year  an  airglow  observing 
station  was  operated  near  Rapid  City,  South  Dakota  (Lat.  ^°35'20", 
Long.  103°l8'50").  A  number  of  observations  were  made  of  stable  auroral 
red  arcs  (hereafter  referred  to  as  SAR-arcs^)  of  the  type  first  described 
by  Barbier  [1958].   Roach  and  Marovich  [i960]  have  described  the  detail- 
ed morphology  of  a  SAR-arc  which  occurred  on  October  22/23,  1958, 
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and  Marovich  and  Roach  [1963]  have  analyzed  the  latitudinal  distribution 
of  23  SAR-arcs  observed  at  Rapid  City. 

A  review  of  the  phenomenon  has  been  given  by  Roach  and  Roach 
[1963];  some  reasonably  well  established  properties  of  these  arcs 
are: 

1.  They  occur  equatorward  of  the  visual  aurorae  during  a  period 
of  increased  magnetic  activity.   From  the  reported  observations, 
they  appear  to  be  oriented  along  parallels  of  geomagnetic  latitude 
and  occur  in  both  hemispheres. 

2.  The  extent  ranges  in  height  from  300-750  km  and  in 
latitude  up  to  ^.5°. 

3.  There  seems  to  be  a  tendency  for  the  intensity  of  a  SAR-arc  to 
remain  quite  stable  for  many  hours  [Rees,  1961],  in  direct  con- 
trast to  the  rapid  fluctuations  characteristic  of  visible 
aurorae . 

k.      They  have  been  observed  at  as  high  a  geomagnetic  latitude  as 
59°  and  as  low  as  kl° .      There  is  some  evidence  that  SAR-arcs  occur 
simultaneously  at  least  halfway  around  the  globe. 

5.  Compared  to  normal  nighttime  conditions,  the  electron 
density  within  the  arc  is  less  and  the  height  of  the  maximum 
density  is  greater. 

6.  The  spectral  composition  consists  predominantly  of  the  63OO- 
636^  A  doublet  of  [oi]. 
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King  and  Roach  [1961]  explored  the  possibility  that  the  excitation 
might  be  due  to  an  enhancement  of  the  rate  of  a  recombination  process  which  ii 
presumed  to  be  the  principal  excitation  mechanism  for  the  nighttime  red 
airglow.   The  reactions  concerned  are: 

0+  +  N2  -  N0+  +  N  (1) 

and 

N0+  +  e  -  N  +  0(1D).  (2) 

The  first  of  these  reactions  is  the  slower  at  altitudes  of  interest  here 
so  that  the  intensity  of  emission  is  proportional  to  the  product  of  0 
and  N    it  was  found,  however,  that  this  would  require  an  increase  in 
Ng  density  by  a  factor  of  nearly  500  to  explain  the  observed  intensity. 

More  recently  several  alternative  possibilities  have  been  dis- 
cussed in  the  literature  by  Megill,  Rees,  and  Droppleman  [1963],  Megill 
and  Carleton  [1964],  Dalgarno  [196W,  and  Cole  [1965].   The  hypotheses 
discussed  in  these  papers  have  in  common  the  requirement  of  excitation 
by  electron  collisions  to  excite  the  0(  D)  but  differ  on  the  method 
by  which  they  propose  to  get  the  energy  into  the  electrons.   They  all 
have  an  implicit  assumption  that  the  arc  is  due  to  an  interaction 
between  the  magnet osphere  and  the  ionosphere. 

The  SAR-arc  phenomenon  is  of  interest  then.,  not  only  in  itself,  but 
as  a  possible  clue  to  mechanisms  operative  during  magnetic  storms.   The 
data  available  in  the  Rapid  City  observations  have  not  been  thoroughly 
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exploited.   As  a  result  the  data  for  23  nights  on  which  the  SAR-arc 
observed  have  been  partially  reduced  and  put  on  magnetic  tape  so  that 
further  investigations  can  be  carried  out  with  a  minimum  of  effort. 
Table  I  gives  a  list  of  the  nights  for  which  data  are  available.  ~ It  is 
the  purpose  of  this  note  to  describe  the  preliminary  operations  which 
have  been  performed  on  these  data  in  placing  them  on  tape  and  to  give 
in  greater  detail  the  behavior  of  a  particularly  good  example  of  the 
arc  on  September  16/17 ,   1958.   These  further  data  involve  the  automatic 
"mapping"  of  the  intensity  in  a  manner  which  may  prove  useful  to  other 
investigators. 

DATA  SEDUCTION 

The  original  data  are  in  the  form  of  strip  chart  records.   Re- 
cords were  obtained  using  a  birefringent  photometer  as  described  by 
Megill  [1959].   The  data  were  obtained  in  five  almucantar  sweeps 
covering  approximately  the  whole  sky  in  five  minutes  for  each  of  the  three 
spectral  lines  observed.   An  entire  sky  survey  was  obtained  every  15 
minutes.   The  field  of  view  of  the  instrument  was  ^°. 

The  data  for  each  degree  in  azimuth  at  zenith  distances  of  80° , 
15° ,   70°,  60°,   k0° ,   and  the  zenith  were  first  placed  on  paper  tape  by 
means  of  a  semi-automatic  reader.   During  this  operation  care  was  taken 
to  exclude  doubtful  data  due  to  the  presence  of  clouds  or  excessive 
Sunlight.   An  approximate  correction  was  made  for  scattering  and 
extinction  in  the  lower  atmosphere  using  the  equation 

I  =R[I-Sc(z)l   ]eTm(z)  (3) 

c  x  '     avg 
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where  I  is  the  record  reading,  z  is  the  zenith  angle,  Sc(z)  is  the 
correction  for  scattering,  t  is  the  extinction  coefficient  and  m(z) 
is  the  air  mass.   (The  air  mass  is  unity  in  the  zenith  at  sea  level). 
The  values  of  Sc(z),  e     'used  are  given  in  Table  II.   The  constant 
E  is  a  calibration  coefficient  dependent  on  the  instrument  and  was  used 
to  convert  the  intensities  to  Rayleighs.   The  average  intensity  is  defined 
by: 

I    =  ±2  [I  -  I  Sc(z)]  eTm  (z)  /V(z)  (10 

avg   n        o  '  x  '  x  ' 

where  n  indicates  the  total  number  of  readings  in  the  survey,  I  is  the 
zenith  intensity,  and  V  is  the  van  Rhijn  function  given  by 

"Vl  -  [E/  (E  +h)]2  sin2  z  (5) 

where  E  is  the  radius  of  the  earth  and  h  is  the  height  of  the  emitting 
layer  above  the  earth's  surface  (taken  as  ^00  km).  What  this  process 
does,  in  effect,  is  to  use  the  zenith  intensity  as  a  first  approximation 
to  the  average  intensity.  A  second  approximation  is  then  obtained  using 
the  deduced  average.   It  should  be  noted  that  the  correction  applies 
strictly  to  a  uniform  intensity  so  that  the  application  of  Eq.  (3)  when 
an  arc  is  in  the  sky  represents  an  approximation  even  if  I    is 
accurately  known.  Using  I    in  performing  the  scattering  calculations 
gives  results  which  are  reliable  to  two  digits  within  the  region  of  the  arc. 
There  appears  to  be  a  tendency  for  this  term  to  overcorrect  in  the  region 
with  no  enhancement  and  the  numbers  are  much  less  reliable.  A  more  detailed 
discussion  of  the  calibration  technique  is  given  in  the  section  on  calibration. 
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"MAPPING"  OF  THE  DATA 
In  order  to  portray  the  data  an  automatic  mapping  program  has  "been 
prepared.   Examples  of  the  output  of  this  program  are  shown  in  the  maps 
included  with  this  note.   This  program  takes  the  data  from  the  tape 
prepared  as  described  above  and  makes  a  first  order  correction  for  the 
van  Rhijn  factor  by  assuming  a  constant  height  layer  at  400  km.   In  order 
to  produce  a  "map"  each  recorded  point  is  assigned  a  coordinate  value 
in  a  polar  coordinate  system.  A  transformation  from  the  polar  coordinate 
values  to  the  map  matrix  is  made  by  assigning  to  each  possible  print 
location  an  intensity  found  by  a  linear  interpolation  between  the 
two  nearest  intensities  on  the  radium  of  the  polar  plot.   All  but  the 
off  scale  values  are  normalized  to  10.25  "the  highest  printed  value.  As 
only  one  digit  per  location  can  be  printed,  values  between  9-75  and 
10.25  are  printed  as  0;  values  between  8.75  and  9.25  are  printed  as  a 
9,   etc.   The  intermediate  values  as  well  as  values  below  .75  are  left 
"blank.  Where  data  were  not  obtained  an  X  is  printed  and  where  the  records 
were  off  scale  a  "+"  is  printed.   The  limit  of  resolution  of  this  system 
varies  from  1  part  in  kO   at  the  high  intensities  to  1  part  in  k   at  the 
lower  intensities.   This  program  allows  the  researcher  to  obtain  a  map 
giving  a  picture  of  the  sky  intensity  without  hand  plotting.   The  night 
of  September  16/YJ ,    1958,  has  been  plotted  in  this  manner  and  included  with 
this  note.   Additional  data  giving  the  time  the  survey  was  taken,  an 
average  intensity  (AVG  INT),  and  the  zenith  intensity  (ZEN  INT)  in 
units  of  Rayleighs  are  printed  at  the  side  of  the  map.   NORM,  when 
multiplied  by  the  digit  of  a  particular  spot  on  the  map,  gives  the 
intensity  in  Rayleighs  for  that  location  (remembering  zero  =10). 
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A  map  relating  the  distance  on  the  earth's  surface  to  the 
corresponding  zenith  angle  is  given  in  Figure  1.   The  height  of  the 
arc  was  assumed  to  "be  ^00  km  throughout  the  analysis.  As  one  would 
expect,  the  height  of  the  SAR-arc  greatly  affects  the  distances  at  the 
large  zenith  angles.   Radial  distances  for  various  heights  have  been 
calculated  as  a  function  of  the  zenith  angle  and  are  shown  in  Figure  2. 

CALIBRATION 

o 

The  absolute  calibration  in  the  6300  A  radiation  rests  on  the 
average  calibration  obtained  for  the  [oi]  (X5577)  radiation  throughout 
the  period  of  operation  of  the  Rapid  City  photometer.   Beginning  with 
the  calibration  for  A-5577  a  factor  was  applied  which  considered  l) 
the  different  photo -multiplier  tube  response,  2)  the  interference  pre- 
filter  characteristics,  and  3)  the  birefringent  filter  characteristics. 
The  A.5577  calibration  in  turn  is  based  on  (a)  intercalibration  with 
the  Fritz  Peak  photometer  by  triangulation  on  the  same  patch  of  night 
sky,  and  (b)  inter -calibration  with  the  Fritz  Peak  instruments  by 
means  of  an  auxiliary  portable  photometer  transported  between  the 
two  stations  on  several  occasions.   The  calibration  of  the  Fritz  Peak 
instruments,  at  the  time,  was  based  on  star  crossings.   The  precision 
(internal  consistency)  of  the  data  is  within  10$.   The  absolute 
calibration  is  probably  not  better  than  a  factor  of  2.   For  a  discussion 
of  accuracy  in  such  photometric  observations  the  reader  is  referred  to 
Roach  and  Smith  [196^-]  and  Smith  and  Alexander  [1963]. 
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SUMMARY 
Preliminary  operations  performed  to  make  the  I.G.Y.  data  on 
Stable  Auroral  Red  Arcs  available  for  quick  usage  are  described.   These 
data  are  available  in  the  form  described  to  researchers  who  wish  to 
use  them.   In  addition,  if  it  should  be  desired,  the  program  which  will 
map  the  data  can  be  made  available.   The  data  are  on  magnetic  tape 
in  binary  format.   The  proper  Fortran  call  and  end  of  record  instructions 
are  given  in  Appendix  I. 


Two  of  the  co-authors,  Davies  and  Rees,  were  supported  in  this 
investigation  by  NSF  Grant  GP620. 
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TABLE  I 
1957  1958 

Nov.  6/7  Jan.  16/17 

Nov.  8/9  Jan.  17/18 

Nov.  9/10  Mar.  Il/l2 

Nov.  lO/ll  Apr.  29/30 

Nov.  25/26  May  27/28 

Nov.  27/28  June  9/l0 

Nov.  30/Dec.  1  June  28/29 

Dec.  lO/ll  Sept.  5/6 

Dec.  20/21  *Sept.  16/17 

Oct.  22/23 
Oct.  23/24 
Oct.  25/26 
Dec.  2/3 
Dec.  17/l8 


*  Data  used  in  this  publication. 
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TABLE  II 


Values  for  6300  A.  [Megill,  1959] 


Zenith 


angle  (z)  Sc(z)  e  *  ' 


0  .059  1.095 

ko  .080  1. 126 

60  .12^  1.200 

TO  .178  1.290 

75  .227  l.kl6 

80  .333  1.665 
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APPENDIX  I 
There  are  two  master  data  tapes  -  one  for  each  year  of  observations. 
The  nights  recorded  on  the  tapes  are  in  the  order  listed  in  Table  I.   By 
using  two  tapes,  the  computer  time  spent  in  searching  for  a  particular 
night  is  decreased  considerably.   It  should  be  noted  that  the  data  have 
been  written  on  the  tapes  in  high  density.  As  they  are  binary  tapes,  no 
format  statement  is  necessary.   The  command  to  recall  this  information  is: 
READ  TAPE  N,  K,  NSTA,  NCLR,  NYR,  MO,  NDAY,  NH,  KM, 
RAVG,  ZEN,  EN,  (jARC(l,J),  1*1,  360),  J  =  1,5^ 


where 


N  -  is  the  designated  tape  unit. 

K  -  a  fixed  point  number  which  is  either  32000  or  zero.  When  a 

zero  occurs,  it  signifies  the  end  of  the  tape.   Unlike  an  end  of 

file,  it  is  possible  for  the  user  to  control  the  tape  at  this  point 

(i.e.,  backspace  and  write,  exit,  etc.). 

NSTA  -  this  is  a  code  used  for  identifying  the  station  which  took 

the  data.   The  number  assigned  the  Eapid  City  Observation  Station 

was  k.      This  number  is  arbitrary  and  may  be  changed  without 

affecting  the  results  of  the  mapping  program. 

o  o 

NCLE  -  As  there  were  three  lines  being  observed  (5577  A>  5893  A, 

o 

and  the  63OO  A),  it  was  necessary  to  identify  the  wavelength  being 

o 

analyzed.   In  this  instance  the  red  line  (6300  A)  was  assigned  the 

o 

fixed  point  number  3.   Although  only  the  6300  A  line  is  presently 
on  the  tape,  the  program  is  set  up  to  include  the  green  line  and 
the  sodium  line.   In  order  to  get  the  correct  printout  for  line 
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identification,  the  number  3  must  be  used  for  the  red  line. 

NYR,  MO,  NBA.Y  -  three  fixed  point  values  which  refer  to  the  year, 

month,  and  the  day  respectively. 

NH,  NM  -  both  of  these  values  are  fixed  point  and  represent  the 

time  of  the  survey  in  hours  and  minutes. 

RAVG  -  the  average  intensity  of  the  sky  after  corrections  for 

extinction  and  scattering  have  been  made. 

ZEN  -  the  zenith  intensity  averaged  over  the  "whole  zenith  sweep. 

EN  -  a  sensitivity  reading  taken  from  the  instrument  at  Rapid  City. 

This  was  for  information  only  and  not  used  in  the  calculations  for 

mapping. 

ARC(l,j)  -  this  array  contains  the  values  of  survey  for  every  degree 

of  azimuth  (I  =  1,360)  for  each  zenith  angle  value  (j  =  1,5).   The 

zenith  angles  correspond  to  the  "J's"  as  follows: 

J  Zen.  Angle 

1  80 

2  75 

3  TO 
k  6o 
5  ko 

Off  scale  values  appear  as  2  x  10  on  the  tape  and  are  printed  as 
+'s.   Unreliable  values  appear  as  -1000  on  the  tape  and  are  printed  as  x's, 
The  original  program  was  written  for  the  IBM  7°90  computer. 
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FIGURE  CAPTIONS 
Figure  1.  A  map  relating  the  distance  on  the  earth's  surface 

to  the  corresponding  zenith  angle  assuming  the  height  of 

emission  is  ij-00  km. 
Figure  2.   Radial  distances  along  the  earth's  surface  for  various 

heights  as  a  function  of  the  zenith  angle. 
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DESCRIPTION  OF  MAPS 

Looking  at  the  map  with  the  legend  in  the  upper  left-hand  corner, 
north  is  at  the  top  of  the  paper  and  east  is  to  the  right.   The  observing 
station  is  at  the  center  of  the  circle.   Zenith  angles  from  the  station 
and  the  corresponding  distances  on  the  earth's  surface  may  be  calculated  as 
in  Figure  1. 

The  arc  is  characterized  by  a  line  of  higher  intensities  from 
«320°  to  %L00°  (0°  at  N  and  90°  at  E).   On  either  side  of  this  line 
the  intensities  are  lower  by  a  substantial  margin.  The  aurora  (first 
quadrant)  is  clearly  discernible  by  the  "+'s"  as  the  intensities  are 
considerably  greater  than  those  of  the  arc. 
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